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Abstract 

This book contains a thorough theoretical consideration of the process of 
proton channeling through carbon nanotubes. We begin with a very brief 
summary of the theoretical and experimental results of studying ion chan- 
neling through nanotubes. Then, the process of ion channeling is described 
briefly. After that, the crystal rainbow effect is introduced. We describe how 
it was discovered, and present the theory of crystal rainbows, as the proper 
theory of ion channeling in crystals and nanotubes. We continue with a de- 
scription of the effect of zero-degree focusing of protons channeled through 
nanotubes. It is shown that the evolution of the angular distribution of chan- 
neled protons with the nanotubc length can be divided in the cycles defined 
by the rainbow effect. Further, we analyze the angular distributions and 
rainbows in proton channeling through nanotubes. This is done using the 
theory of crystal rainbows. The angular distributions are generated by the 
computer simulation method, and the corresponding rainbow patterns are 
obtained in a precise analysis of the mapping of the impact parameter plane 
to the transmission angle plane. We demonstrate that the rainbows enable 
the full explanation of the angular distributions. The analysis demonstrates 
that the angular distributions contain the information on the transverse lat- 
tice structure of the bundle. In this study we investigate the rainbows in 
transmission of protons of the kinetic energy of 1 GeV through a straight 
very short bundle of (10, 10) single-wall carbon nanotubc, a bent short bun- 
dle of (10, 10) nanotubes, and the straight long (11, 9) nanotubes. 

We also investigate how the effect of dynamic polarization of the carbon 
atoms valence electrons influences the angular and spatial distributions of 
protons transmitted through short nanotubes in vacuum and embedded in 
dielectric media. It is demonstrated that this effect can induce the additional 
rainbow maxima in the angular distributions. We have established that the 
changing of the spatial distribution with the proton kinetic energy, based 
on the changing of the dynamic polarization effect, may be used to probe 
the atoms or molecules intercalated in the nanotubes. In these study we 
investigate the rainbows in transmission of 0.223-2.49 MeV protons through 
the short (11, 9) nanotubes in vacuum and embedded in Si02, AI2O3 and 
Ni. Besides, we report on the donut effect in transmission of 0.223 MeV 
protons through a short (11, 9) nanotube, which occurs when the initial ion 
velocity vector is not parallel to the nantoube axis. In addition, we explore 
the channeling star effect in 1 GeV proton channeling through bundles of 
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nanotubes, which appears when the proton beam divergence angle is larger 
than the critical angle for channeling. 



1 Introduction 



Carbon nanotubes were discovered by lijima in 1991 (jlijimal . Il99ll ). One can 
describe them as the sheets of carbon atoms rolled up into the cylinders with 
the atoms lying at the hexagonal crystal lattice sites. The diameters of nan- 
otubes are of the order of a nanometer and their lengths can be more then a 
hundred micrometers. Nanotubes can be the single-wall and multi-wall ones, 
depending on the number of cylinder s they include. Th ey have remarkable 



geometrical and physical properties ( ISaito et al.l . 1200 ll ). As a result, nan- 



otubes have begun to play an important role in the field of nanotechnologies. 
For example, they are expected to beco me the basic elements for creating 



nanoelectronic devices (lYao et al.l . Il999l ) 



The structure of a carbon nanotube is described by vector Ch = mai + 
na2 = {m,n), where ai and 02 are the unit vectors of the hexa gonal crystal 
lattice, and the integers satisfying inequalities < \n\ < m (ISaito et al. 



2OO1I ). This vector is called the chiral vector of the nanotube. Its magnitude 



is Ch = {jn^ + ran + n'^)2ah, where ah is the magnitude of ai and 02. The 
angle of Ch relative to ai, which is called the chiral angle of the nanotube, is 
given by expression cosB/i = {m + n/2)/(m^ + mn + 71^)2 and inequalities 
< |6h| < 7r/6. This is shown in Fig. [T] The sheet of carbon atoms is rolled 
up into the cylinder in such a way to transform the segment corresponding to 
Ch into a circle. Hence, the nanotube diameter equals Ch/'n'- The nanotube 
is achiral or chiral. If it is achiral, the nanotube consists of the atomic strings 
parallel to its axis. In particular, if n = 0, when B/^ = 0, or m = n, when Qh 
= 7r/6, the nanotube is achiral. In accordance with the shape of its transverse 
cross-section, in the former case the nanotube is called the zigzag one, and 
in the latter case the armchair one. If it is chiral, the nanotube consists of 
the atomic strings that s piral around its a xis. This is illustrated in Fig. [2l 

fcooih . 



Saito et al 



which is taken from Ref. 

Soon after the d i scover y of carbon nanotubes, Klimov and Letokhov 
(IKlimov fc Letokhovl . Il996l ) predicted that they could be used to channel 
positively charged particles. After that, a number of theoretical groups have 



studied ion channeling through rianotubes (jGevorgian et al.l. ll998l:IZhevago fc Glebov 



19981 : iBirvukov fc Bellucci 120021 : iGreenenko fc Shulgal . l2003l : IZhevago fc Glebov 
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2005[lMoura fc Amarai ^2005':'Bork a et al.l. 120051: iNeskovic et al.l . j2005l:IPetrovic et al. . 



2005a 



b; Borka et al 



2006aiib. ,20071: iMatvuk hin fc Frolenko^. 120071: iMiskovid 



20071: iMoura fc Amarail2007l : lBorka et al.l.l2008a.b 



2008al Jbl: IZheng et al.l . l2008aU bl: iPetrovic et al.l . l2009l ). The main objective of 



Matvukhinl . 120081 : iPetrovic et al. 



most of those studies was to investigate t he possibility of guid i ng ion beams 
with nanotubes. Biryukov and Bellucci (IBiryukov &: Belluccil . |2002| ) looked 
for the nanotube diameter opti mal for channeling high energy ion beams. 
Krasheninnikov and Nordlund (IKrasheninnikov fc Nordlund] , |2005| ) studied 
the channeling of low energy Ar"*" ions through achiral and chiral nanotubes. 
Petrovic et al. (IPetrovic et al.l . l2005al ) demonstrated that the rainbow effect 
could play an important role in ion channeling through nanotubes. 




Figure 1: The hexagonal crystal lattice: xq and yo are the axes of the local 
reference frame, and ai and 02 the unit vectors of the lattice. When carbon 
atoms are placed at the lattice sites, a sheet of carbon atoms is obtained. 
When rolled up into a cylinder, the sheet of carbon atoms becomes a car- 
bon nanotube: Ch is the chiral vector of the nanotube, Qh the chiral angle 
of the nanotube, and T the translational vector of the nanotube, which is 
perpendicular to Ch and extends from its origin (point O) to the first lattice 
point. Rectangle defined by Ch and T is a unit cell of the nanotube. In the 
displayed case Ch = (6,8) and T = (4,-5). 
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Figure 2: The classification of carbon nanotubes: (a) an armcliair nanotube, 
(b) a zigzag nanotube, and (c) a cliiral nanotubes. 



Tlie experimental studying of ion channeling through carbon nanotubes 
is in the initial phase. The most challenging task in such experiments still is 
to solve the problems of ordering, straightening and holding nanotubes. The 
first experiment al data on i on ch anneling through nanotubes were reported 
by Zhu et al. (IZhu et al.l . |2005| ). They were obtained with He"^ ions of 



the kinetic energy of 2 MeV and an array of the well-ordered multi-wall 
nanotubes grown in a porous anodic aluminum oxide (AI2O3) membrane. 
The authors measured and compared the yields of ions transmitted through 
the bare AI2O3 sample and the AI2O3 sample including the nanotubes. The 
first experiment with electrons and nanotubes was performed by Chai et 
al. (IChai et al.l . 120071 ). They used the 300 keV electrons and studied their 
transport through the multi-wall nanotubes of the length in the range of 
0.7-3.0 fim embedded in a carbon fiber c oating. The misalignm ent of the 



nanotubes was below 1°. Berdinsky et al. (IBerdinsky et al.l . |2008[ ) succeeded 



in growing the single-wall nanotubes in the etched ion tracks in a silicon 
oxide (8102) layer on a silicon substrate. This result has opened a possibility 
for conducting precise measurements of ion channeling through nanotubes in 
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a wide ion energy range. 

This chapter is devoted to our theoretical studies of proton channehng 
through carbon nanotubes. We shall first describe briefly the process of ion 
channeling. Then, the crystal rainbow effect will be introduced, and the 
theory of crystal rainbows will be presented. We shall continue with a de- 
scription of the effect of zero-degree focusing of protons channeled through 
nanotubes. Further, the angular distributions and rainbows in proton chan- 
neling through nanotubes will be analyzed. We shall focus on the rainbows 
in transmission of 1 GeV protons through a straight very short bundle of (10, 
10) single-wall carbon nanotube, a bent short bundle of (10, 10) nanotubes, 
and the straight long (11, 9) nanotubes. Also, the influence of the effect of 
dynamic polarization of the carbon atoms valence electrons on the angular 
and spatial distributions of protons transmitted through short nanotubes in 
vacuum and embedded in dielectric media will be investigated. We shall con- 
centrate on the rainbows in transmission of 0.223-2.49 MeV protons through 
the short (11, 9) nanotubes in vacuum and embedded in Si02, AI2O3 and 
Ni. Besides, we shall report on the donut effect in transmission of 0.223 MeV 
protons through a short (11, 9) nanotube, which occurs when the initial ion 
velocity vector is not parallel to the nantoube axis. In addition, we shall 
explore the channeling star effect in 1 GeV proton channeling through bun- 
dles of nanotubes, which appears when the proton beam divergence angle is 
larger than the critical angle for channeling. 



2 Ion channeling 



An axial crystal channel is a part of a crystal in between its neighboring 
atomic strings being parallel to one of its crystallographic axes. The process 
of ion motion through the axial crystal channel, in which the angle of its 
velocity vec tor relative to the channel axis remains small, is called ax ial ion 
channeling (IRobinson and Peru . Il963l : iLindhardl . Il965t iGemmelll . Il974j ) . This 
process is explained by the ion repulsions from the atomic strings deflning 
the channel, which are the results of the series of its correlated collisions 
with the atoms of the strings. Analogously, a planar crystal channel is a 
part of the crystal in between its neighboring atomic planes being parallel 
to one of its crystallographic planes. The process of ion motion through the 
planar crystal channe l is called plana r ion channeling (IRobinson and Oenl . 
19631 : ILindhardl . Il965l : iGemmell Il974l ). 
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Theoretical studies of ion channeling through crystals have bee n going on 
along two major lines. The first line was founded by Lindhard ( iLindhardl . 



19651 ) and the second one by Barrett (IBarrettl . Il97ll ). The Lindhard's ap- 



proach was analytical and it was developed by the methods of statistical 
mechanics. He included the continuum approximation, i.e., neglected the 
longitudinal correlations between the positions of the atoms of one atomic 
string, did not take into account the transverse correlations between the 
positions of the atomic strings, and included the assumption of statistical 
equilibrium in the transverse position plane. The Barrett's approach was via 
the ion-atom scattering theory and it was numerical. He developed a very 
realistic computer code for the three-dimensional following of ion trajectories 
through crystal channels, in which the three approximations used by Lind- 
hard were avoided. The Barrett's approach is more complex, but a number 
of calculations have shown that it is much more accurate than the Lindhard's 
approach. 

The results that will be presented here have been obtained by a realis- 
tic computer code based on the numerical solution of the ion equations of 
motion through crystal channels. The continuum approximation is included. 
Consequently, one has to solve only the ion equations of motion in the trans- 
verse position plane. The transverse correlations between the positions of 
the atomic strings are taken into account. 

Let us take that the z axis of the reference frame coincides with the 
channel axis and that its origin lies in the entrance plane of the crystal. The 
X and y axes of the reference frame are the vertical and horizontal axes, 
respectively. It is ass umed that the interaction of the ion and cryst al can 
be treated classically (ILindhardl . Il965l : iGemmelll . Il974j : iBarrettl . Il97ll ). The 
interaction of the ion and a crystal's atom is described by an approximation of 
the Thomas-Fermi interaction potential or any other appropriate interaction 
potential - V{r'), where r' is the distance between the ion and crystal's atom. 
As a result, we obtain the ion-crystal continuum interaction potential. 



J 

U{x,y) = J2Uj{x,y), (1) 

where 
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+00 

U,{x,y) = ^ I V{p^ + z')Uz, (2) 
—00 

is the continuum interaction potential of the ion and the j*^ atomic string, 



Pi 



[{x - X,)' + {y- 



1/2 



(3) 



is the distance between the ion and this atomic string, and x and y are the 
transverse components of the proton position vector; J is the number of 
atomic strings and d the distance between the atoms of an atomic string. 

The thermal vibrations of the crystal's atoms can be taken into ac- 
count. This is done by averaging U{x, y) over the tran sverse displacera e nts of 



the crystal's at oms from their equilibrium positions (lAppleton et al.l . 11967 : 



Neskovid . Il986l ). As a result, one has to substitute Uj{x,y) in Eq. ([2]) with 



th/ 



a 



X, y) = Uj{x, y) + — [d^^Uj{x, y) + dyyUj{x, y)] 



(4) 



where ath is the one-dimensional thermal vibration amplitude of the crystal's 
atoms. Thus, U^^{x,y) is obtained instead of U{x,y). 

We neglect the nuclear ion energy loss, i.e., the energy loss resulting from 
its collisions with the crystal's nuclei. But, the electronic ion energy loss and 
dispersion of its transmission angle, which are caused by its collisions with 
the crystal's electrons, can be included. For the specific ion energy loss we 
use expression 



dE 

dz 



2^4 



AnZte 



-Hf, 



In 



2me7^t'(-2)^ 

hUf. 



(5) 



where Zi is the atomic number of the ion, e the elementary charge, rrie the 
electron mass, v{z) the magnitude of the ion velocity vector, (3 = v{z)/c, 



c the speed of light, 7 



He 



{d^x + dyy)U"'{x,y)/An the 



average of the density of the crystal's electron gas along the z axis, h the 
reduced Planck constant, and Ue = {ine^rie/me)^ the ang ular freque r icy of 



the ion induced osci llations of the crystal's electron gas (iGemmelll . 11974 : 



Petrovic et al.l . |2002| ). For the specific change of the dispersion of the ion 



transmission angle we use expression 
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where m = 777,07 relati v istic ion mass and mo the ion rest mass 



(iGemmelll . Il974j : iPetrovic et aLl . l2002l ). This dispersion is a measure of the 
uncertainty of the ion transmission angle. The dispersions of the vertical 
and horizontal components of the ion transmission angle, B^; and 0^, are 
dnl^ = dVLl/2 and dVtl^ = dVLl/2, respectively. 

The ion equations of motion in the transverse position plane are 



777 



d'^x 
dF 



and 



th. 



(7) 



777 



dt^ 



x,y), 



where t denot es the time. They are solved by the Runge-Kutta method of the 
fourth order (jPress et al.l . Il993[ l. In the calculation the time step is chosen 
in such a way to obtain the fixed longitudinal distance step with the desired 
accuracy. 

At the end of each time step of the calculation the magnitude of the ion 
velocity vector is corrected using the value of the ion energy loss during the 
step. The vertical and horizontal components of the ion transmission angle 
at the end of each time step are chosen within the Gaussian distribution 
functions defined by 6^. and dQ^^ and by 9^ and dQ^y, respectively. The 
(relativistic) ion mass at the end of each time step is calculated using the 
corrected value of the magnitude of the ion velocity vector. 

The initial direction of the ion velocity vector is defined by the angles 
its vertical and horizontal components make with the x and y axes, Q^q and 
Qyo, respectively. However, if this direction is uncertain, with the dispersions 
of Qxo and G^o being Q"^^ = and Q'^y = ^^^2, respectively, where fid 

is the divergence angle of the initial ion beam, the vertical and horizontal 
components of the initial ion velocity vector are chosen via the Gaussian dis- 
tribution functions defined by Q^o and fi^^ and by G^o and fl'^y, respectively. 



3 Crystal rainbows 

Let us now examine the process of ion transmission through an axial crystal 
channel. The scheme of the process is given in Figure El It is well-known that 
the motion of an ion close to the channel axis can be treated as an oscillatory 
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motion around the axis. If the ion impact parameter vector and the angle 
between its initial velocity vector and the channel axis arc fixed, the number 
of oscillations the ion makes before leaving the crystal depends on its initial 
kinetic energy. If this energy is sufficiently high for the ion to make less 
than about a quarter of an oscillation around the channel axis, its trajectory 
can be approximated by a straight line. If this is true for the majority of 
ions, one says that the crystal is very thin. If, however, the majority of ions 
make between about a quarter of an oscillation and about one oscillation, 
the crystal is thin. The crystal is thick if the majority of ions make more 
than about one oscillation. Finally, if the majority of ions make much more 
than one oscillation, the crystal is very thick. 




Figure 3: A scheme of the process of ion transmission through an axial crystal 
channel. The z axis of the reference frame coincides with the channel axis 
and its origin lies in the entrance plane of the crystal. The x and y axes are 
the vertical and horizontal axes, respectively. In this example the channel is 
defined by four atomic strings of the crystal, the initial ion velocity vector 
lies in the xz plane and it is parallel to the z axis. The vertical components 
of the ion impact parameter vector and transmission angle are designated by 
X and Qx, respectively. The horizontal components of these variables, y and 
Qy, respectively, are equal to zero. 
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3.1 Discovery of crystal rainbows 

Neskovic f lNeskovicl . [l986h showed theoretically that in the ion transmission 
through an axial channel of a very thin crystal the rainbow occurred, i.e., the 
ion differential transmission cross section could be singular. His approach was 
via the ion-molecule scattering theory. In the created model the continuum 
approximation was included implicitly, the correlations between the positions 
of the atomic strings of the crystal were included, and the assumption of sta- 
tistical equilibrium in the transverse position plane was avoid ed. The effect 



was o bserved experimentally for the first time by Krause et al. (iKrause et al 



1986 ). It was named the crystal rainbow effect ( Neskovic et al. . 19931 ). 



The first measurement of crystal rainbows was perfor med with 7 M e V pro 



tons and the <100> and <110> Si very thin crystals (IKrause et al.l . Il986l ). 



The obtained angular distrib utions of chann eled protons were reproduced 
using the model of Neskovic ( Neskovic . 19861) wi t h the Lindhard's proton- 
crystal continuum interaction potential (ILindhardl . Il965l ) . Krause et al. have 
also performed the high resolution measurements of the angular distribu- 
tions of 2-9 MeV protons and 6-30 MeV C^"*", C^"* " and C^"*" ions transmitted 
through the <100> Si thin crystal ( IKrause et al.l . Il994l ). The interpretation 



of the obtained results was done using the computer code of Barrett (iBarrett 



19711 ) with the Moliere's ion-crystal continuum interaction potential. Those 



measurements demonstrated clearly the possibilities of using crystal rainbows 
for exploring the properties of thin crystals, e.g., for precise local measure- 
ment of their thickness, with the accuracy of a few atomic layers, and for 
precise differential me asurement of the electron d ensity and ion energy loss 
within their channels (INeskovic fc Petrovig . 120031 ) . 



As it has been said in part §1 of this chapter, Petrovic et al. ( iPetrovic et al. 



2005al ) showed theoretically that the rainbow effect can also appear with car- 
bon nanotubes. This will be discussed in detail in part §5 of this chapter. 



3.2 Theory of crystal rainbows 



Krau se et al. (IKrause et al.l . Il986l ) and Miletic et al. ( iMiletic et al.l . Il996l . 



19971 ) investigated the periodicity of evolution of the angular distribution of 
channeled ions with the reduced crystal thickness, which is 



A 



f{q,m)L 



(9) 
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where q = Zie, m and vq are the ion charge, mass and magnitude of the 
initial ion velocity vector, respectively, L is the crystal thickness, and f{q, m) 
is the average frequency of the ion motion close to the channel axis. They 
found that the evolution of the angular distribution of channeled ions could 
be divided into cycles. The first cycle lasts for A between and 0.5, the 
second cycle for A between 0.5 and 1, and so on. These cycles were named 
the rainbow cycles. We note that the crystal is very thin if A is smaller than 
~0.25, it is thin if A is between ~0.25 and ~1, it is thick if A is larger than 
~1, and it is very thick if A is much larger than ~1. For A = (n — l)/2 + 1/4, 
where n = 1, 2, ... is the order of the rainbow cycle, a large part of the ion 
beam is focused around the origin in the transverse position plane. This is an 
effec t of spatial focusin g ; whic h is called the eff e ct sup erfocusing of channeled 



ions 
= 1, 



Demkov fc Meverl fl2004[ ): iNeskovic et al.l (120091). For A = n/2, with n 



2, a large part of the ion beam is parallel to the channel axis, i.e., 
it is focused around the origin in the transmission angle plane. This is an 
effect of angular f ocusing, which is called th e effect of zero-degree focusing 
of channeled ions ( iMiletic et al.l . Il996l . 119971 ) . 

In order to describe the ion transmission through axial channels of crys- 
tals that are not necessarily very t hin, Petrovic et al. ( jPetrovic et al.l . |2000|) 
generalized the model of Neskovic ( iNeskovid . Il986l ) and formulated the the- 
ory of crystal rainbows. Let us describe this theory briefly. We take that the 
z axis of the reference frame coincides with the channel axis and its origin 
lies in the entrance plane of the crystal. The x and y axes of the reference 
frame are the vertical and horizontal axes, respectively. This is shown in Fig. 
|3l It is ass umed tha t the interaction of the ion and crystal can be treated 
classically ( iLindhardl . Il965l ). T he interact i on of the ion and a crystal's atom 
i s described by the Lindhard's ( ILindhardl . Il965l ) or Moliere's approximation 
(iMolierd . Il947l ) of the Thomas-Fermi interaction potential, or any other ap- 
propriate interaction potential. We apply the continuum approximation, i.e., 
the atomic string s are treated as if they are continual, rather than discrete 
(ILindhardl . Il965l ). As a result, we obtain the Lindhard's or Moliere's ion- 
crystal continuum interaction potential, or another appropriate continuum 
interaction potential. The thermal vibrations of the crystal 's atoms can be 



taken into account as it is described in part 2 of this chapter ( lAppleton et al. 



19671 ). We neglect the nuclear ion energy loss, i.e., the energy loss resulting 



from its collisions with the crystal's nuclei. But, the electronic ion energy loss 
and dispersion of its transmission angle, which are caused by its collisions 
with the crystal's electrons, can be included as it is described in part 2 of 



14 



D. Borka, S. Petrovic and N. Neskovic 



this chapter (iGemmelll . ll974J : iPetrovic et al.l . 120021 ) 



The transverse components of the initial ion position vector are equal to 
the components of its impact parameter vector, xq and yo, while the trans- 
verse components of the initial ion velocity vector are v^o and VyQ. In order 
to obtain the components of the final ion position and velocity vectors, one 
must solve the ion equations of motion. Since we apply the continuum ap- 
proximation and the electronic ion energy loss during the whole transmission 
process remains small, we assume that the longitudinal ion motion is uniform 
and solve only the ion equations of motion in the transverse position plane. 
In order to compensate the inaccuracies resulting from this assumption, the 
magnitude of the ion velocity vector at the end of each step of the calculation 
is corrected using the value of the ion energy loss during the step. We also 
correct the direction of the ion velocity vector at the end of each step of the 
calculation using the value of the change of the dispersion of its transmission 
angle during the step. This is how the ion energy loss and change of the 
dispersion of its transmission angle are included in the calculations. The rel- 
ativistic effect is taken into account by changing the ion mass at the end of 
each step of the calculation using the corrected value of the magnitude of the 
ion velocity vector. The obtained transverse components of the final ion po- 
sition vector are x{xo, yo, A) and y{xo, yo, A), and the transverse components 
of the final ion velocity vector are Vx{xo, yo, A) and Vy{xo, yo, A). Since the ion 
transmission angle during the whole transmission process remains small, the 
components of the final ion transmission angle are Qx{xo,yo, A.) = v^/v and 
Qy{xo, yo, A) = Vy/v, where v = v{xo, yo. A) is the magnitude of the final ion 
velocity vector. The spatial and angular distributions of transmitted ions, 
in the exit plane of the crystal, are generated by the computer simulation 
method. The spatial and angular distributions of ions in the entrance plane 
of the crystal are determined by the spatial and angular characteristics of 
the chosen initial ion beam. If the ion, during its motion along the channel, 
enters any of the cylindrical regions around the atomic strings of the radius 
equal to the screening radius of a crystal's atom, it is excluded from the 
calculations. 

Since the ion transmission angle during the whole transmission process 
remains small, the ion differential transmission cross section is 



a{xo, yo, A) = — — (10) 

\j0{xo,yo,A)\ 



where 
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Je(xo, Vo, A) = ^^o^x^yo©?/ - dyoQ^d^aQy (11) 

is the Jacobian of the components of its transmission angle, which equals 
the ratio of the infinitesimal surfaces in the transmission angle plane and 
impact parameter plane. It describes the mapping of the impact parameter 
plane to the transmission angle plane. Hence, equation jQ{xo,yo, A) = 
gives the angular rainbow lines in the impact parameter plane, along which 
the mapping is singular. The images of these lines determined by functions 
Qx{xo, yo, A) and Qy{xo, yo, A) are the rainbow lines in the transmission angle 
plane. 

On the other hand, the mapping of the impact parameter plane to the 
transverse position plane is described by the Jacobian of the transverse com- 
ponents of the ion position vector. This Jacobian, which equals the ratio of 
the infinitesimal surfaces in the transverse position plane and impact param- 
eter plane, reads 



Jp{xo, yo, A) = d^^xdy^y - dy^^xd^^y 



(12) 



Thus, equation Jp{xo, yo, A) = gives the spatial rainbow lines in the impact 
parameter plane, along which the mapping is singular. The images of these 
lines determined by functions x{xo,yo,A.) and y{xo,yo,A.) are the rainbow 
lines in the transverse position plane. 

The rainbow lines in the transverse position plane and transmission angle 
plane separate the bright and dark regions in these planes. Their shapes 
are classified by catas trophe theory ( jThoml . Il975t iNeskovic fc Perovig . Il987 : 



Neskovic et al.l . 119931 ). Hence, one says that the ion beam dynamics in the 



channel has the catastrophic character. 



4 Zero-degree focusing of protons channeled 
through bundles of nanotubes 

The first calculations of the effect of zero-degree focusing of channeled ions, 
which is defined in part ^3 .2 of this chapter, were performed by Miletic et al. 
(IMiletic et al.l . Il996l . 119971 ). In the former study the projectiles were 25 MeV 
C^^ ions and the target was the <100> Si crystal of the thickness between 
and 2.7 /im, while in the latter study the projectiles were Ne^°+ ions of the 
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kinetic energies between 1 and 37 MeV while the target was the 1.0 fim thick 
<100> Si crystaL The effect has not yet been observed experimentally. 
We shall describe here the effect of zero-degree focusing o f 1 GeV protons 



channeled through the (10, 1 0) single-wall carbon nanotubes (INeskovic et al 



20051 : IPetrovic et aP . l2005bh . Such a nantoube is achiral, i.e., its atomic 



strings are parallel to its axis. It is assumed that the nanotubes form a bun- 
dle whose transverse cross-section can be descirbed via a (two-dimens ional) 



hexagonal or rhombic superlattice with one nanotube per primitive cell (iThess et al. 



We choose the z axis of the reference frame to coincide with the bundle 
axis and its origin to lie in the entrance plane of the bundle. The arrange- 
ment of the nanotubes is such that their axes intersect the x and y axes of 
the reference frame, w hich are the vertical and horizontal axes, respectively 



(IPetrovic et al.l . l2005al ). We take into account the contributions of the nan- 



otubes lying on the two nearest rhombic coordination lines, relative to the 
center of the primitive cell of the (rhombic) superlattice. 

The calculations are performed using the theory of crystal rainbows, 
which is described in part §3.2 of this chapter. The interaction of th e proton 
and a nanotube atom is described by the Moliere's approximation (IMoliere , 



19471 ) of the Thomas-Fermi interaction potential, which reads 



V{r') 



ZiZ2e 



[0.35exp(-0.3r7a) + 0.55 exp(-1.2r7a) + 0.10 exp(-6.0r7a)] 

(13) 

1 and Z2 = Q are the atomic numbers of the proton and nanotube 



is the distance between 



where Zi 

atom, respectively, e is the elementary charge 

1 

the proton and nanotube atom, b = 0.3/a, a = [97r^/(128Z2)] ^ Oq is the 
screening radius of the nanotube atom, and Cq is the Bohr radius. It has been 
proven that this expression provides excellent agreemen t with experimental 



results in the field of ion channeling (IKrause et al.l . Il994( ). The application of 



the continuum approximation gives the Moliere's proton-bundle continuum 
interaction potential, which reads 



(14) 



=1 i=i 



where 
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2Z Z 

Uij{x,y) = y [0.35iro(0.3p,,/a) + ().bbK^{l.2pi^/a) + 0.10iro(6.0p,,/a)] 

(15) 

is the Moliere's continuum interaction potential of the proton and the jth 
atomic string of the i^^ nanotube within the bundle, 



Pij = [{x - Xijf + {y- Vij)^] 



(16) 



is the distance between the proton and this atomic string, and x and y are the 
transverse components of the proton position vector; / = 16 is the number 
of nanotubes within the bundle, J = 40 is the number of atomic strings of 
a nanotube, d = 0.24 nm is the distance between the atoms of an atomic 
string, and Kq denotes the modified Bessel function of the second kind and 
Qth Qj-fjgr xhis means that the number of atomic strings within the bundle 
is / X J = 640. The thermal vibrations of the nanotube atoms are taken into 
account. This is done by substituting Uij{x,y) in Eq. (IT^ with 



(x, y) = U,,{x, y) + -^ [d^^U,,{x, y) + dyyUij{x, y)] 



(17) 



where ath = 5.3 p m is the one-dimen sional thermal vibration amplitude of the 



nanotube atoms ( Hone et al. . 2000l ). Thus, U^'^{x,y) is obtained instead of 



U{x,y). The electronic proton energy loss and dispersion of its transmission 
angle, caused by its collisions with the nanotube electrons, are neglected. 
The nanotube radius is 0.67 nm ( ISaito et al.l . l200ll ) and the distance be tween 



the axes of two neighboring nanotubes is 1.70 nm ( iThess et al.l . Il996l ). The 



components of the proton impact parameter vector are chosen uniformly 
within the primitive cell of the (rhombic) superlattice. The initial proton 
velocity vectors are all taken to be parallel to the bundle axis. 

The nanotube walls define two separate regions in the transverse position 



plane: inside the nanotubes and in between them (iPetrovic et al.l . I2005al ). 



Accordingly, the bundle contains two types of channels: the circular one, 
whose center coincides with the center of the region inside each nanotube, 
and the triangular one, whose center coincides with the center of the region 
in between each three neighboring nanotubes. It is clear that in this case 
one has to define two average frequencies of the proton motion close to the 
channel axis, corresponding to the protons moving close to the centers of 
the circular and triangular channels. The two frequencies can be determined 
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from the second order terms of the Taylor expansions of the proton-bundle 
continuum interaction potential in the vicinities of the centers of the two 
types of channels. Consequently, there are two reduced bundle lengths, Ai 
and A2, corresponding to the proton motions close to the axes of the circular 
and triangular channels, respectively. 

Figure IH^a) shows the dependence of the zero-degree yield of protons 
propagating along the circular channels of the bundle on the bundle length, 
L, in the range of 0-200 /im; the Ai axis is shown too. For the region in 
the transmission angle plane in the vicinity of its origin we take the region 
in which the proton transmission angle, G = (©^ + ©y)^, where G^; and 
Qy are its vertical and horizontal components, respectively, is smaller than 
0.0109 mrad. The initial number of protons is 174 976. The six maxima of 
the dependence correspond to the ends of the first six rainbow cycles, where 
the proton beam channele d through the circular channels is quasi-parallel 



(IMiletic et al.l . Il996l 119971 ). The positions of these maxima should be close 



to Ai = 0.5, 1, 1.5, .... However, this is not true. For example, the position 
of the first maximum of the dependence is Ai = 0.35. The large deviations of 
the values of these position from their expected values are explained by the 
strong anharmonicity of the proton-bundle continuum interaction potentia l 



in the vicinities of the centers of the circular channels (IMiletic et al.l . 119971 ). 



The dependence of the zero-degree yield of protons propagating along the 
triangular channels of the bundle on L in the range of 0-200 /im is given in 
Fig. m^b); the A2 axis is shown too. The initial number of protons is 136 658. 
The dependence has 27 maxima, and they correspond to the ends of the first 
27 rainbow cycles, where t he proton beam channele d through the triangular 



channels is quasi-parallel ( IMiletic et al.l . Il996l Il997l ). The positions of these 
maxima are close to A2 = 0.5, 1, 1.5, .... For example, the position of the 
first maximum of the dependence is A2 = 0.49. The small deviations of the 
values of these position from their expected values are attributed to the weak 
anharmonicity of the proton-bundle continuum int eraction potential in the 



vicinities of the centers of the triangular channels (IMiletic et al.l . Il997l ). It 



should be noted that the total zero-degree yield of channeled protons, i.e., 
of protons channeled through both the circular and triangular channels, as 
a function of L is the sum of the zero-degree yields shown in Fig. IH^a) and 
Fig. Hb). 

We have also investigated the effect of zero-degree focusing o f 1 GeV pro- 
tons channeled through the (10, 0) single-wall carbon nanotubes ( Borka et al.l . 
20081 ) . The study has been performed in the same way as in the case of (10, 
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Figure 4: (a) The zero-degree yield of 1 GeV protons channeled through 
the circular channels of the bundle of (10, 10) single-wall carbon nanotubes 
as a function of the bundle length, L. (b) The zero-degree yield of 1 GeV 
protons channeled through the triangular channels of the bundle of (10, 10) 
single-wall carbon nanotubes as a function of the bundle length, L. 



10) nanotubes. The number of nanotubes within the bundle is / = 16, the 
number of atomic strings of a nanotube is J = 20, the distance between the 
atoms of an atomi c string is d = 0.21 nm. The nanotube radius is 0.39 nm 
(ISaito et all . 1200 ll ). an d the distance betw een the axes of two neighboring 
nanotubes is 1.10 nm (jZhang et al.l . |2003| ). In this case the dependence of 



the zero-degree yield of protons channeled through the circular channels of 
the bundle on L in the range of 0-200 fim contains 26 maxima. The posi- 
tions of these maxima are close to Ai = 0.5, 1, 1.5, .... The average distance 
between the maxima is AAi = 0.49. This means that the anharmonicity 
of the proton-bundle continuum interacti on potential in the v icinities of the 



centers of the circular channels is weak ( iMiletic et al.l . 119971 ). In this case 



the dependence of the zero-degree yield of protons channeled through the 
triangular channels of the bundle on L in the range of 0-200 /xm contains 37 
maxima. The positions of these maxima are very close to A2 = 0.5, 1, 1.5, 
.... The average distance between the maxima is AA2 = 0.50. Thus, one 
can say that the anharmonicity of the proton-bundle continuum interaction 
potential in the vicinities of the centers of the triangular channels is very 
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weak flMiletic et allll997h . 



The comparison of the zero-degree yields of protons channeled through 
the circular and triangular channels of the bundles as functions of L in the 
cases of (10, 10) and (10, 0) nanotubes shows that: 



- (i) when the circular channels are considered, the average frequency of 
the proton motion close to the channel axis is much lower in the former 
than in the latter case, 

- (ii) when the triangular channels are considered, the average frequency 
of the proton motion close to the channel axis is lower in the former 
than in the latter case, and 

- (iii) when the circular and triangular channels are compared to each 
other, the ratio of the average frequencies of the proton motions close 
to the channel axes in the circular and triangular channels is much 
smaller in the former than in the latter case. 



Thus, one can conclude that the measurements of the effect of zero-degree 
focusing of channeled protons can give the information on the transverse 
lattice structures of the nanotubes and of the bundles. 



5 Angular distributions and rainbows in pro- 
ton channeling through nanotubes 

We shall present in this part of the chapter the angular distributions and 
rainbows in channeling 1 GeV protons t hrough a straight very short bundle 
of (10, 10) single-wall carbon n anotubes flPetroyic et al.l . l2005al ). a bent short 



bundle of (10, 10) nanotubes (Neskovic et al. , 20051 ). and the straight long 



(11, 9) nanotubes f lPetrovic et al.l . l2008bf ) 



5.1 Rainbows with a straight very short bundle of na- 
notubes 



The system we investigate here is a 1 GeV proton moving through a bundle 
of (10, 10) single-wall carbon nanotu bes that is desc r ibed in part 4 of this 
chapter. The bundle length is 1 fim ( iPetrovic et al.l . l2005al ). The reduced 



Angular distributions and rainbows. 



21 



bundle lengths corresponding to the proton channeling inside the nanotubes 
and in between them, i.e., through the circular and triangular channels of 
the bundle, are Ai = 0.015 and A2 = 0.070, respectively. These values tell 
us that in both cases the majority of protons make before leaving the bundle 
less than a quarter of an oscillation around the channel axis (Ai, A2 < 0.25). 
Therefore, one can say that the bundle we investigate is very short. 




-1 

y (nm) 



Figure 5: The rainbow lines in the impact parameter plane for 1 GeV protons 
transmitted through the 1 /im long (10, 10) single-wall carbon nanotubes. 

The calculations are performed using the theory of crystal rainbows, 
which is described in part §3.2 of this chapter. The interaction of the proton 
and bundle is described by Eqs. (fT4l) -( fT6l) . We take into account the effect 
of thermal vibrations of the nanotube atoms, and this is done by Eq. ( |T71) . 
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Figure 6: The rainbow lines in the transmission angle plane corresponding to 
the rainbow lines in the impact parameter plane (shown in Fig. 1) lying (a) 
inside each nanotube, and (b) in between each four neighboring nanotubes. 



However, the electronic proton energy loss and dispersion of its transmission 
angle are neglected. This is justified by the fact that the bundle is very short. 
If the bundle were long (Ai,A2 > 1), these effect wou ld cause the smeariri g 



out of the angular distributions of channeled protons ( jPetrovic et al.l . l2002l ). 



Figure [5] shows the rainbow lines in the impact parameter plane in the 
case under consideration. One can see that inside each nanotube, i.e., inside 
each circular channel of the bundle, there is one (closed) rainbow line, while 
in between each three neighboring nanotubes, i.e., inside each triangular 
channel of the bundle, there are one larger and four smaller (closed) rainbow 
lines. 

The rainbow line in the transmission angle plane that is the image of 
the rainbow line in the impact parameter plane lying inside each nanotube 
is shown in Figure [6]^a). It consists of 20 connected cusped triangular lines 
lying along the lines ip = tan~^(Gy/6a;) = 2{n + l)7r/20, n = 0-19, which 
correspond to the parts of the rainbow line in the impact parameter plane 
in front of the 20 pairs of atomic strings defining the nanotube (see Fig. 
[5]). Points 1 and 2 are the intersection points of the rainbow line in the 
transmission angle plane with the line Qy = 0. Points 1 are the apices of the 
cusps and points 2 are the intersections of the parts of the rainbow line. The 
corresponding points in the impact parameter plane are also designated by 
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Figure 7: (a) The angular distribution of 1 GeV protons transmitted through 
the 1 fim long bundle of (10, 10) single-wall carbon nanotubes. The areas in 
which the yields of transmitted protons are larger than 0.13, 0.26, 0.39, 1.3, 
2.6, 3.9, 13, 26, 39 % of the maximal yield are designated by the increasing 
tones of gray color, (b) The corresponding yield of transmitted protons along 
the line 6^ = 0. 



1 and 2. 

Figure Mjo) shows the rainbow lines in the transmission angle plane that 
are the images of the rainbow lines in the impact parameter plane lying 
in between each four neighboring nanotubes. The analysis shows that the 
rainbow pattern consists of two cusped equilateral triangular rainbow lines 
in the central region of the transmission angle plane with the cusps lying 
along the lines cp = 2nn/3 and = {2n + l)vr/3, n = 0-2, each connected 
with three pairs of cusped triangular rainbow lines lying along the same lines, 
and eight cusped triangular rainbow lines lying in between the six pairs of 
triangular lines. The two equilateral triangular lines each connected with the 
three pairs of triangular lines are the images of the two larger rainbow lines 
while the eight triangular lines are the images of the eight smaller rainbow 
lines in the impact parameter plane (see Fig. [5]). Points 1', 2', 3' and 4' are 
the intersection points of the rainbow lines in the scattering angle plane with 
line Qy = 0. Points 1' are the apices of the cusps, and points 2', 3' and 4' are 
the intersections of the parts of the larger rainbow lines. The corresponding 
points in the impact parameter plane are also designated by 1', 2', 3' and 4'. 

Figure [7](a) shows the angular distribution of transmitted protons. The 
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number of transmitted protons is 2 142 538. The areas in which the yields of 
transmitted protons are larger then 0.13, 0.26 and 0.39 %, 1.3, 2.6 and 3.9 %, 
and 13, 26 and 39 % of the maximal yield are designated by the increasing 
tones of gray color. At the very low level of the yield, corresponding to the 
boundary yields of 0.13, 0.26 and 0.39 % of the maximal yield, there are 
20 triangular forms in the peripheral region of the transmission angle plane, 
with the maxima lying on the lines = 2(n + l)7r/20, n = 0-19. Further, 
at the low level of the yield, corresponding to the boundary yields of 1.3, 2.6 
and 3.9 % of the maximal yield, there is a hexagonal structure in the central 
region of the transmission angle plane, with the maxima lying on the lines 
(p = mr/S, n = 0-5. Finally, at the high level of the yield, corresponding 
to the boundary yields of 13, 26 and 39 % of the maximal yield, there is 
a pronounced maximum at the origin of the transmission angle plane. The 
analysis shows that the first part of the angular distribution is generated by 
the protons with the impact parameters close to the atomic strings defining 
the nanotubes - the 20 triangular forms correspond to the 20 pairs of atomic 
strings defining the nanotubes. The second part of the angular distribution 
is generated by the protons with the impact parameters in between the nan- 
otubes but not close to the centers of the triangular channels. The third part 
of the angular distribution is generated to a larger extent by the protons 
with the impact parameters close to the centers of the circular channels and 
to a smaller extent by the protons with the impact parameters close to the 
centers of the triangular channels. It should be noted that most of the pro- 
tons that generate the third part of the angular distribution interact with the 
nanotubes very weakly - they move through the space inside the nanotubes 
virtually as through a drift space. Thus, we can say that the angular dis- 
tribution contains the information on the transverse lattice structure of the 
bundle. Its first part (the peripheral region of the transmission angle plane) 
provides the information on the individual nanotubes while its second part 
(the central region of the transmission angle plane) provides the information 
on the way they are connected to each other. 

The comparison of Figs. El^a) and MJo) with Fig. [Tl^a) shows clearly 
that the shape of the rainbow pattern determines the shape of the angular 
distribution of transmitted protons. Also, each maximum of the angular 
distribution, except the maximum lying at the origin of the transmission 
angle plane, can be attributed to one of the above mentioned characteristic 
rainbow points in the transmission angle plane. Thus, one can conclude that 
the rainbow pattern enables the full explanation of the angular distribution. 
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Figure H^b) gives the low and very low levels of the yield of transmitted 
protons along line Gy = 0. The arrows indicate the above mentioned char- 
acteristic rainbow points in the transmission angle plane. It is evident that 
the two maxima at the low level of the yield can be explained by points 1' 
and 4' in the transmission angle plane, and the two shoulders at the very low 
level of the yield by points 1, 2, 2' and 3'. This means that the characteristic 
rainbow points in the transmission angle plane can provide the information 
on the continuum potential of the bundle at the corresponding points in the 
impact parameter plane, and, hence, on the average electron density in the 
bundle at these points (see Fig. [5]). The two maxima at the low level of the 
yield can be used to measure the average electron density at points 1' and 4', 
lying in between the nanotubes, while the two shoulders at the very low level 
of the yield can be used to measure the average electron density at points 
1 and 2, lying in the nanotube, and at points 2' and 3', lying in between 
the nanotubes. The obtained data can help one compare various theoretical 
approaches and determine the electron structure of the bundle. 



We have also performed the analyses of the angular distributions of 1 
GeV protons transmitted through the 1 /im long bundles of (10, 0) and (5, 
5) single-wall carbon nanotubes. These nanotubes are achiral too. In the 
former case each nanotube consists of 20 atomic strings and in the latter 
case of 10 pairs of atomic strings. The angular distribution in the (10, 0) 
case is similar to the one in the (10, 10) case while the angular distribution in 
the (5, 5) case is very different from the one in the (10, 10) case. In both cases 
it is easy to establish the correspondence between the parts of the angular 
distributions and the transverse lattice structures of the bundles. It should 
be also noted that t he evolution of the angular distribution in the (10, 0) case 
(IBorka et al.l . 120051 ) with the proton en ergy or the bundle len gth is different 
from the evolution in the (10, 10) case flPetrovic et"aD . l2005bl ). enabling one 
to distinguish between the two types of bundles. Thus, the obtained results 
can lead to a new method of characterization of achiral carbon nanotubes, 
based on the rainbow effect. This method would be complementary to the 
existing method of characterization of nanotubes by electrons impinging on 
them transverse ly rather than longi tudinally, which is based on the diffraction 



effect (see, e.g., iLucas et al.l fl2002[ )) 
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5.2 Rainbows with a bent short bundle of nanotubes 

Again, the system we explore is a 1 GeV proton moving through a bundle 
of (10, 10) single-wall carbon nanotubes that is described in part §4 of this 
chapter. However, now, the bundle is bent along the y axis. The bundle 
l ength is 7 /im and it s bending angle is varied in the range of 0-1.5 mrad 
(INeskovic et al.l . l2005l ). In this case the critical angle for channeling is \l/c = 



0.314 mrad. If the angle of the proton velocity vector relative to the channel 
axis is above "^C: the proton is dechanneled. The reduced bundle lengths 
corresponding to the proton channeling through the circular and triangular 
channels of the bundle are Ai = 0.103 and A2 = 0.488, respectively. Thus, in 
the circular channels the majority of protons make before leaving the bundle 
less than a quarter of an oscillation around the channel axis (Ai < 0.25), and 
in the triangular channel they make close to half an oscillation (A2 ~ 0.5). 
Hence, for the protons in the circular channels the bundle is very short while 
for those in the triangular channels it is short (Ai < 1). 

The calculations are performed using the theory of crystal rainbows, 
which is described in part §3.2 of this chapter. The interaction of the proton 
and bundle is described by Eqs. f|T^ - f|T6l) . We take into account the effect 
of thermal vibrations of the nanotube atoms, and this is done by Eq. f|T71) . 
However, the electronic proton energy loss and dispersion of its transmission 
angle are neglected. This is justified by the above given values of Ai and A2. 

The bundle length has been chosen to correspond to the first maximum 
of the zero-degree yield of protons propagating along the triangular channels 
of the bundle as a function of the bundle length, which is given in Fig. 
111(b). Thus, the part of the proton beam transmitted through the triangular 
channels is quasi-parallel, and we can investigate in an easier manner the 
behavior of the proton beam transmitted through the circular channels. 

Figure 5.4(a) shows the angular distribution of protons transmitted through 
the straight bundle, i.e., for its bending angle a = 0. The sizes of a bin along 
the Qx and Qy axes are equal to 0.00667 mrad while the initial number of 
protons is 866 976. The angular distribution contains (i) a pronounced max- 
imum at the origin of the transmission angle plane, (ii) six non-pronounced 
maxima lying along the lines defined by = ta.n~^{Qy/Qx) = nir/S, n = 
0-5, and (iii) a non-pronounced circular part. The rainbow pattern in the 
transmission angle plane obtained for a = is given in Fig. (Ht^b). It contains 
(i) two six-cusp lines very close to the origin of the transmission angle plane, 
corresponding to two complex lines in the impact parameter plane close to 
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Figure 8: (a) The angular distribution of 1 GeV protons transmitted through 
the straight bundle of (10, 10) single-wall carbon nanotubes of the length of 
L = 7 fim. (b) The corresponding rainbow pattern in the transmission angle 
plane. 



the centers of the two triangular channels making each rhombic channel of 
the bundle, (ii) six points lying along the lines defined by = ravr/S, n = 
0-5, corresponding to two times three points in the impact parameter plane 
within the two triangular channels making each rhombic channel, and (iii) a 
circular line, corresponding to a circular line in the impact parameter plane 
within each circular channel. The comparison of the results given in Figs. 
(Ht^a) and 5.4(b) shows that the angular distribution can be explained as fol- 
lows: (i) the pronounced maximum at the origin of the transmission angle 
plane is generated by the protons with the impact parameter vectors close 
to the centers of each circular channel and each triangular channel, (ii) the 
six non-pronounced maxima by the protons with the impact parameter vec- 
tors close to the six rainbow points (within each rhombic channel), and (iii) 
the non-pronounced circular part by the protons with the impact parameter 
vectors close to the circular rainbow line (within each circular channel). 

Figure [9](a) shows the angular distribution of protons transmitted through 
the bundle obtained for a = 0.2 mrad. The transmission angle plane is 
parallel to the exit plane of the bundle. The sizes of a bin along the 
and Qy axes and the initial number of protons are the same as in the case 
of a = 0. One can see easily that the angular distribution lies in the region 
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Gy > -0.2 mrad, demonstrating that the proton beam is bent by the bundle 
effectively. The centrifugal force made the average bending angle of the 
proton beam smaller than 0.2 mrad. The angular distribution contains (i) 
a part having the shape of an acorn with two maxima at points (±0.096 
mrad, -0.054 mrad), and (ii) eight additional maxima close and in between 
the two maxima. The rainbow pattern in the transmission angle plane for 
a = 0.2 mrad is given in Fig. MJo). It contains (i) an acorn-like line with 
two joining points of its branches, corresponding to a complex line in the 
impact parameter plane within each circular channel of the bundle, and (ii) 
eight groups of points, corresponding to two times four groups of points in 
the impact parameter plane within the two triangular channels making each 
rhombic channel. The comparison of the results given in Figs. M^a) and 
Mjo) shows that the angular distribution can be explained as follows: (i) the 
acorn-like part with the two maxima is generated by the protons with the 
impact parameter vectors close to the acorn-like rainbow line with the two 
joining points (within each circular channel), and (ii) the eight additional 
maxima by the protons with the impact parameter vectors close to the eight 
groups of rainbow points (within each rhombic channel). 

Figure [TD] gives the dependences of the yields of protons transmitted 
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Figure 9: (a) The angular distribution of 1 GeV protons transmitted through 
the bent bundle of (10, 10) single-wall carbon nanotubes of the length of L 
= 7 fim for the bending angle a = 0.2 mrad. (b) The correspoding rainbow 
pattern in the transmission angle plane. 
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Figure 10: The yields of protons transmitted through the circular channels, 
the triangular channels, and both the circular and triangular channels of the 
bent bundle of (10, 10) single- wall carbon nanotubes of the lenght of L = 7 
/xm as functions of the bending angle, a. 



through the circular channel and the triangular channels of the bundle on a 
in the range of 0-1.5 mrad. The initial numbers of protons are 77 866 for 
the former dependence and 60 770 for the latter dependence. The inflection 
points of the dependences appear at a = 0.35 and 0.45 mrad, respectively. 
These values ought to be compared to the value of which is 0.314 mrad. 
However, the former yield decreases with a after the inflection point consid- 
erably slower than the latter yield. They fall below 10 % of the yields for 
the straight bundle at a = 1.45 and 0.75 mrad, respectively. This figure also 
gives the total yield of protons transmitted through the bundle, both through 
the circular and triangular channels, as a function of a. Its inflection point 
appears at a = 0.45 mrad, and it falls below 10 % of the yield for the straight 
bundle at o; = 1.25 mrad. 
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5.3 Rainbows with long nanotubes 



Let us now consid er a 1 GeV proton mov ing; through a long (11, 9) single- wall 
carbon nanotube (jPetrovic et all l2008bl ). The nanotube length, L, is varied 
from 10 to 500 fim. Such a nanotube is chiral, i.e. its atomic strings spiral 
around its axis. The z axis of the reference frame coincides with the nanotube 
axis and its origin lies in the entrance plane of the nanotube. The x and y 
axes of the reference frame are the vertical and horizontal axes, respectively. 

In the calculations we employ the theory of crystal rainbows, which is 
decribed in part §3.2 of the chapter. The interaction of th e proton and a 
nanotube atom is described by the Moliere's approximation (iMolierd . 119471 ) 
of the Thomas- Fermi interaction potential, which is given by Eq. (fT3l) . Since 
the nanotube is chiral, the interaction potential of the proton and nanotube 
is obtained by the azimuthal averaging of the Moliere's proton-nanotube 
continuum interaction potential, which is given by Eqs. ffT^ - ffTB]) . with the 
thermal vibrations of the nanotube atoms taken in to account, via Eq . (d 
in these equations I = 1, J 
interaction potential reads 



40, and ath = 5.3 pm (iThess et al.l . ll996l ). This 
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wher e acc = 0.14 nm is the bond length o f the carbon atoms (jSaito et al. 
20011 ) ■ R = 0.69 nm is the nanotube radius (ISaito et al.l . 1200 ll ). (ctj) = (0.35, 

(x^ + 



y 

and 



0.55, 0.10) and (/3j) = (0.1, 1.2, 6.0) are the fitting parameters, p 
'^^2^ X and y are the transverse components of the proton position vector. 
In deno t es th e modified Bessel function of the first kind and O*'* order 
Artru et al.l (120051 ). It is evident that U^^{x,y) is cylindrically symmetric. 

The electronic proton energy loss and dispersion of its transmission an- 
gle, caused by its collisions with the nanotube electrons, are not taken into 
account. The components of the proton impact parameter vector are chosen 
randomly within the circle around the origin of radius R — a. The initial 
proton velocity vectors are all taken to be parallel to the nanotube axis. 

Since U^^{x,y) is cylindrically symmetric, the problem we consider is 
in fact one-dimensional. This means that the rainbow lines in the impact 
parameter plane and transmission angle plane, if they exist, will show up as 
the circles. Consequently, it is sufficient to analyze the mapping of the x 
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axis in the impact parameter plane to the G axis in the transmission angle 
plane, i.e., the Qx{x) deflection function. The abscissas and ordinates of the 
extrema of this deflection function determine the radii of the rainbow lines 
in the impact parameter plane and transmission angle plane, respectively. 

Figures [T]T a)-(d) show the angular distributions of channeled protons 
along the B axis for L = 10, 50, 100 and 500 fim, respectively. The size 
of a bin along the G axis is 0.866 /xrad and the initial number of protons is 
16 656 140. These values enabled us to attain a high resolution of the angular 
distributions in a reasonable computational time. The angular distributions 
obtained for L = 10, 50 and 100 /im contain a central maximum and a number 
of symmetric pairs of maxima characterized by a sharp decrease of the proton 
yield on the large angle side; the numbers of symmetric pairs of sharp maxima 
are one, eight and 15, respectively. The angular distribution obtained for L = 
500 yum contains a central maximum and a large number of symmetric pairs 
of sharp maxima close to each other. The analysis shows that in the region 
where IG^^I < 0.2 mrad one can identify 55 pairs of sharp maxima. Figure 
[TlT d) also shows that in the region where 0.041 mrad < |Ga;| < 0.058 mrad 
there are flve pairs of sharp maxima; they are designated by 11-15. In the 
region in which [G^, ! > 0.2 mrad the resolution of the angular distribution is 
not sufficiently high to distinguish easily between the adjacent pairs of sharp 
maxima. 

The Ga;(x) deflection functions obtained for L = 10, 50, 100 and 500 m 
are given in Figs. [T^ a)-(d). respectively. The deflection functions for L = 
10, 50 and 100 /xm contain one, eight and 15 symmetric pairs of extrema, 
respectively; each pair includes a minimum and a maximum. The analysis 
shows that the deflection function obtained for L = 500 fim contains 76 
symmetric pairs of extrema. Figure [T^ d) also gives flve pairs of extrema, 
designated by 11-15, that correspond to the flve sharp maxima shown in Fig. 
[TlT d). Each symmetric pair of extrema deflnes a circular rainbow line in 
the impact parameter plane and a circular rainbow line in the transmission 
angle plane. One can observe that each of the deflection functions for L = 
50, 100 and 500 /zm has two envelopes, one of them connecting its extrema 
designated by odd numbers and the other connecting its extrema designated 
by even numbers, and that it oscillates between these envelopes. 

The comparison of Figs. [TT](a)-(c) and Figs. [T2r a)-(c) shows that the ab- 
scissas of the pairs of sharp maxima of each angular distribution of channeled 
protons along the G^^ axis coincide with the ordinates of the corresponding 
pairs of extrema of the accompanying Qx{x) deflection function. One can 
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Figure 11: The angular distributions along the G^; axis of 1 GeV protons 
channeled through the (11, 9) SWCNTs of the length of (a) 10 /im, (b) 50 
/im, (c) 100 yUm and (d) 500 /im. 



draw the same conclusion for the pairs of sharp maxima of the angular dis- 
tribution given in Fig. [TlT d) and the corresponding pairs of extrema of the 
deflection function given in Fig. [T2r d). in spite of the fact that in the region 
where 10^:1 > 0.2 mrad it is not possible to identify the additional 21 pairs of 
sharp maxima of the angular distribution. Hence, it is clear that the sharp 
maxima of the angular distributions can be attributed to the rainbow effect. 
The bright and dark sides of the rainbows are the small and large angle sides 
of the corresponding sharp maxima. 

Figures [T2t^a)-(d) show clearly that the number of rainbows generated by 
the transmitted protons increases and the average distance between them 
decreases as L increases. These two dependences, for L between 50 and 
500 fim, are given in Fig. [121 The numbers and positions of the rain- 
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Figure 12: The deflection functions G^; of the x axis in the impact parameter 
plane for 1 GeV protons channeled through the (11, 9) SWCNTs of the length 
of (a) 10 /im, (b) 50 /im, (c) 100 /im and (d) 500 fim. 



bows are determined from the accompanying Qx{x) deflection functions. 
The analysis shows that the former dependence can be approximated ex- 
cellently by fitting function , where oi = 0.15 is the fitting parameter. 
The latter dependence can be approximated excellently by fitting function 
/2 = a2exp(— L/62) + 03 exp(— L/63), where 02 = 0.072 mrad, 62 = 42.6 /im, 
as = 0.015 mrad and 63 = 340.5 /xm are the fitting parameters. These fitting 
functions are also shown in Fig. [131 When the nanotube becomes sufficiently 
long for the average distance between the rainbows to become smaller than 
the resolution of the angular distribution, one cannot distinguish between 
the adjacent rainbows. This means that the rainbows disappear and the an- 
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Figure 13: The dependences of the number of circular rainbows in the angular 
distributions of 1 GeV protons channeled through the (11, 9) single- wall 
carbon nanotubes - N, open circles, and the average distance between them 
- S, closed circles, on the nanotube length - L. The full lines represent the 
fitting curves. 



gular distribution becomes a bell-shaped one. The analysis shows that, in 
parallel, the spatial distribution of channeled protons in the exit plane of the 
nanotube also becomes a bell-shaped one. Hence, one can say that, when 
the nanotube becomes sufficiently long, the angular and spatial distributions 
equilibrate, and, as one wo uld expe c t, thi s does not happen in accordance 
with the ergodic hypothesis iBarrettI (Il973l ). We call this route to equilibra- 
tion, which is characterized by the linear increase of the number of rainbows 
and the exponential decrease of the distance between them as L increases, 
the rainbow route to equilibration. 

Recently, a theoretical analysis of the channeling of 1 GeV pr otons through 



the lon g (10, 10) single- wall carbon nanotubes has been reported lPetrovic et al. 



([20091). The angular distributions of transmitted protons were generated for 
L between 10 and 100 /xm. The obtained results show that when L < 30 
Hm the transverse lattice structure of the nanotube can be deduced from the 
angular distribution. When L > 40 fim the angular distribution contains the 
concentric circular ridges whose number increases linearly and the distance 
between them decreases exponentially as L increases, in a similar way as in 
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the above described case of (11, 9) nanotubes. Consequently, these ridges 
can be attributed to the rainbow effect. 

The above presented resuhs demonstrate that one can employ the an- 
gular distributions of protons channeled through long carbon nanotubes for 
their characterization. Namely, it seems that each type of nanotube for each 
value of L has a characteristic pattern of concentric circular ridges, and that 
this can be checked experimentally. However, an easier way to check this 
experimentally would be to vary the proton kinetic energy rather than L. 
This method of characterization of nanotubes would be complementary to 
the method proposed in part §5.1 of this chapter. 



6 Dynamic polarization effect in proton chan- 
neling through short nanotubes 

In this part of the chapter we shall investigate how the angular and spatial 
distributions of protons channeled through a short (11, 9) single- wall carbon 
nanotube in vacuum or embedded in a dielectric medium is influenced by 
the effect of dynana i c polar ization of the nanotube atoms valence electrons 
fiBorka et al.l . l2006aL[2008bl Jah. The magnaitude of the proton velocity vector. 



V, is varaied between 3 and 8 a.u., corresponding to the proton kinetic energy 
between 0.223 and 1.59 MeV, respectively. It is expected that the dynamic 
polarization effect, which is induced by the proton, is pronounced in this 
range of its kinetic energy. The nanotube length, L, is varied between 0.1 and 
0.8 fim. We shall also explore the influence of the angle of the initial proton 
velocity vector relative to the nanotube axis on the angular distribution of 
protons channele d through a short ( 11, 9) nanotube in vacuum for v = 3 a.u. 



and L = 0.2 /im flBorka et al.l . l2010f l. 

The z axis of the reference frame coincides with the nanotube axis and 
its origin lies in the entrance plane of the nanotube. The x and y axes of 
the reference frame are the vertical and horizontal axes, respectively. The 
calculations are performed using the theory of crystal rainbows, which is 
described in part §3.2 of the chapter. The thermal vibrations of the nanotube 
atoms are not taken into account. The electronic proton energy loss and 
dispersion of its transmission angle, caused by its collisions with the nanotube 
electrons, are neglected too. The components of the proton impact parameter 
vector are chosen randomly within the circle around the origin of radius R—a, 



36 



D. Borka, S. Petrovic and N. Neskovic 



where R is the nanotube radius, a = [97r^/(128Z2)]^ cto the screening radius 
of the nanotube atom, and ao the Bohr radius. The initial proton velocity 
vectors are all taken to be parallel to the nanotube axis. 



6.1 Rainbows with short nanotubes in vacuum 

The system under consideration here is a proton moving through a short (11, 
9) single-wall carbon nanotube in vacuum for v = 3 a.u. and L = 0.1-0.3 fim 
(iBorka et al.l . l2006al ) : the corresponding proton kinetic energy is 0.223 MeV. 
This nanotube is chiral. 

The interaction of the proton and a nanotube ato r n is d escribed by 
the Doyle- Turner interaction potential (IDoyle fc Turnerl . Il968l ). Since the 
nanotube is chiral, the static part of the interaction potential of the pro- 
ton and nanotube is obtained by the azimuthal averagi ng of the Doyle- 



Turner proton-naiiotube contin uum interaction potential (iLindhardl . 11965 



Zhevago &: Glebovl . Il998l . |2000| ): the number of atomic strings of the nan- 
otube is J = 40. This interaction potential, in atomic units, is 



Us{x,y) 



'iVoO'cc i=i 



(19) 



where (7^) = (0.115, 0.188, 0.072, 0.020) and (5^) = (0.547, 0.989, 1.982, 
5.656) are the fitting parameters (in a.u.), p = (x^ + 1/^)2, x and y are 
the transverse components of the proton position vector, and Jq denotes the 
modified Bessel function of the first kind and O*'' order. It is clear that 
Us{x,y) is repulsive and cylindrically symmetric. 

The dynamic polarization effect is treated by a two-dimensional hydrody- 
namic model of the nanotube atoms valence electro ns based on a je l lium-li ke 
description of the io n cores on the nanotube walls ( Mowbray et al. . 2004al lbl: 



Borka et al. 



2006ah . This model includes the axial and azimuthal averaging 
as in the case of Us{x,y). As a result, the dynamic part of the interaction 
potential of the proton and nanotube is 



00 

uAx,y)=^f:p[ irf(^i?)"^^±i;Mj^.(,,),fc, (20) 



where 
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u;f{k) = + /2/^2) [47inoRIi{kR)Ki{kR) + v^,] (21) 

is the square of the angular frequency of the proton induced oscillations of the 
nanotube electron gas of angular mode / and longitudinal wave number k, no 
= 0.428 the ground state density of the nanotube electron gas, Vs = (vrno)2 
is the velocity of propagation of the density perturbations of the nanotube 
electron gas, < p < R, Ii and Ki denote the modified Bessel function of 
the first and second kinds and order, respectively, and P designates that 
only the principal part of the integral is taken into account. One can see that 
Ud{x,y), which is also called the image interaction potential of the proton 
and nanotube, is attractive and cyllindrically symmetric. Thus, the total 
interaction potential of the proton and nanotube is 

U{x,y) = Us{x,y) + Ud{x,y). (22) 
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Figure 14: The total interaction potential of the proton and nanotube, 
U{x,y), along the x axis in the interval {—R,R) for v = 3 a.u - solid line. 
The static and dynamic parts of the interaction potential are also shown - 
dotted and dashed lines, respectively. 



Figure [T3] shows the total interaction potential of the proton a nd nanotube , 
U{x,y), along the x axis in interval {—R,R) for f = 3 a.u. ( iBorka et al.l . 
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2006al : iPetrovic et al.l . l2008al ). The static and dynamic parts of the inter- 
action potential are also shown in the figure. One can see that due to the 
inclusion of the image interaction potential the total interaction potential 
of the proton and nanotube has two infiection points. They appear for x 
= ±8.1 a.u. If the nanotube we consider is short, the corresponding Qx{x) 
defiection function will have a pair of symmetric extrema appearing for the 
values of x close to ±8.1 a.u. 
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Figure 15: (a) The angular distributions along the B^^ axis of protons trans- 
mitted through a (11, 9) single-wall carbon nanotubes for v = 3 a.u. and 
L = 0.1 fim with and without the dynamic polarization effect included - 
solid and dashed lines, respectively, (b) The corresponding 6^(a;) deflections 
functions. 



The angular distribution of transmitted proto ns along the 0-^ axi s for v 
= 3 a.u. and L = 0.1 fim is given in Fig. [TST a) (IBorka et al.l . l2006al ). This 
flgure also shows the corresponding angular distribution when the image in- 
teraction potential of the proton and nanotube is not included. The size of 
a bin along the G^. axis is equal to 0.066 mrad and the initial number of 
protons is 3 141 929. One can see that a central maximum and a pair of 
symmetric peripheral maxima appear in the angular distribution when the 
image interaction potential is included; the peripheral maxima are labeled 
by 1. Figure ITHb) shows the Qx{x) deflection functions corresponding to the 
angular distributions given in Fig. [TST a). It is evident that the deflection 
function without the image interaction potential taken into account has no 
extrema. On the other hand, the deflection function with the image interac- 
tion potential taken into account exhibits a pair of symmetric extrema; they 
are also labeled by 1. The analysis shows that the two peripheral maxima 
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of the angular distribution are connected to the two extrema of the deflec- 
tion function. Therefore, we conclude that they are the rainbow maxima, and 
that they occur as a consequence of the inclusion of the dynamic polarization 
effect. Besides, one should note that the extrema of the deflection function 
appear for x = ±8.0 a.u. Since these positions are close to the positions of 
the inflection points of the total interaction potential, one should conclude 
that the nanotube in question is short. 

We have established that for v = 3 a.u. and L = 0.2 /xm the angular 
distribution of transmitted protons along the B^; axis contains three pairs 
of rainbow maxima, and that for v = 3 a.u. and L = 0.3 nm the an gular 



distribution includes five pairs of rainbow maxima (IBorka et al.l |2006aJ) . All 



these maxima are attributed to the image interaction potential of the proton 
and nanotube. 



6.2 Rainbows with short nanotube embedded in di- 
electric media 

Experimentally, the main problem that prevent the observation of ion chan- 
neling through carbon nanotubes is the problem of ordering, straightening 
and holding nanotubes. The most promising method enabling us to solve 
this problem is b ased on the growing of nanotubes inside the holes in a di- 



electric medium (jZhu et al.l . l2005l : iBerdinsky et al.l . |2008|). Besides, in many 



applications of nanotubes it is desirable to have them embedded in SiO^ 



(|Tsetseris fc Pantelidesl . 120061 ). or clamped by a Ni shield (iGuerret-Piecourt et al. 
19941 : iMiskovicl . l2007f ). 



When a nanotube is embedded in a dielectric medium, the static part of 
its interaction potential with the proton is the same as in the case when it 
is in vacuum, which is given by Eq. (19). The dynamic part of the interac- 
tion potential is obtained by a hydrodynamic model of the nanotube atoms 
valence electrons based on the jellium-like description of the ion cores on the 
nanotube walls, as in the case when the nanotube is in vacuum, but extended 
to include the polar ization of the dielectric boundary (IMowbray et al.l . 12007 : 
Borka et al.l . l2008bl ). It should be noted that this model has also been ap- 
plied successfully in the r elated area of ion interactio n with dielectric suported 



sheets of carbon atoms (iRadovic et al. 
Figure ITOT a 



20071 . l2009l ). 



shows the angular distribution along the G^. axis of protons 
transmitted through a (11, 9) single-wall carbon nanotube embedded in Si02 
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for V = 5 a.u. and L = 0.5 /im (iBorka et al.l . l2008bf l: the corresponding pro- 
ton kinetic energy is 0.619 MeV. This figure also shows the corresponding 
angular distribution when the nanotube is in vacuum. The size of a bin 
along the Qx axis is equal to 0.0213 mrad and the initial number of protons 
is 3 141 929. One can see a central maximum and only one pair of sym- 
metric peripheral maxima, labeled by 1^;, when the nanotube is embedded in 
Si02, in comparison with a central maximum and three pairs of symmetric 
peripheral maxima, labeled by 1, 2' and 2" , when the nanotube is in vacuum. 
Figure [T6]^b) gives the Qxix) defiection functions corresponding to the an- 
gular distributions given in Fig. irGT a). We have found that each peripheral 
maximum of each angular distribution is connected to an extremum of the 
corresponding defiection function. Hence, all the peripheral maxima are the 
rainbow maxima. 

Let us know analyze the spatial distributions along the x axis in the exit 
plane of the nanotube of protons transmitted through a (11, 9) single- wall 
carbon na notube embedded in Si02 and in vacuum for v = 5 a.u. and L 
= 0.5 fim ( iBorka et al.l . l2008al ). It is given in Fig. HEJ^c). The size of a bin 
along the 0^; axis is equal to 0.3 a.u. and the initial number of protons is 
3 141 929. In the former case the spatial distribution constains a central 
maximum and three pairs of symmetric peripheral maxima, labeled by 1^, 
2rf and 3^, and in the latter case the spatial distribution includes a central 
maximum and four pairs of symmetric peripheral maxima, labeled by Ij, 2j, 
3j and 4j. The maxima designated by 2^, 3d, 3j and 4j are very weak. Figure 
fTBT d) gives the mappings of the xq axis in the entrance plane of the nanotube 
to the X axis in its exit plane corresponding to the spatial distributions given 
in Fig. flGT c). It is evident that when the nanotube is embedded in Si02, 
the mapping has six symmetric extrema, labeled by 1^^, 2^ and 3^. When the 
nanotube is in vacuum, the mapping has eight symmetric extrema, labeled 
by li, 2j, 3i and 4j. Extrema 2^, 3^, 3j and 4j are very sharp and lie near the 
nanotube wall. In each case the ordinate of the extremum of the mapping 
coincides with the abscissa of the corresponding peripheral maximum of the 
spatial distribution. Therefore, these maxima are attributed to the rainbow 
effect. 

Figure [TTT a) gives the angular distribution along the axis of protons 
channeled through a (11, 9) single-wall carbon nanotube embe dded in Si02 



AI2O.S and Ni and in vacuum for f = 8 a.u. and L = 0.8 fim ( IBorka et al 



2008b|); the corresponding proton kinetic energy is 1.59 MeV. The size of 



a bin along the axis is equal to 0.0213 mrad and the initial number of 
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protons is 3 141 929. The angular distributions obtained with the embed- 
ded nanotube contain one central maximum and two symmetric peripheral 
maxima, labeled by 1^. It is clear that these positions do not depend much 
on the type of medium surrounding the nanotube. This can be explained 
by the fact that in these cases the depend ences of the image force acting 



on the proton on x are close to each other (IBorka et al.l . |2008b|) . When the 



nanotube is in vacuum, the angular distribution contains a central maxi- 
mum and no symmetric peripheral maxima. Figure [T71(b) shows the Qx{x) 
deflection functions corresponding to the angular distributions shown in Fig. 
[TTT a). It is clear that in the cases of embedded nanotubes the deflection 
functions contain two pairs of symmetric extrema. The analysis shows that 
the two peripheral maxima appearing in the angular distributions obtained 
with the embedded nanotube are connected to the former pair of extrema of 
the corresponding deflection functions, also labeled by 1^. This means that 
these maxima are due to the rainbow effect. The angular distributions ob- 
tained with the embedded nanotube do not include any maxima that would 
be connected to the latter pair of extrema of the corresponding deflection 
functions, lying close to the nanotube wall. It is also clear that in the case 
of nanotube in vacuum the deflection function includes only one pair of sym- 
metric extrema, lying close to the nanotube wall. However, as it has been 
already said, the angular distribution obtained with the nanotube in vaccum 
does not have any maxima that would be connected to the pair of extrema 
of the corresponding deflection function. 

Now, we are going to consider the spatial distributions along the x axis in 
the exit plane of the nanotube of protons channeled through a (11, 9) single- 
wall carbon nanotube embed ded in 8102, ^ ^120.^ and Ni and in vacuum for 
V = 8 a.u. and L = 0.8 fim (IBorka et al.l . l2008al ). They are shown in Fig. 
[TTT c). The size of a bin along the G^^ axis is equal to 0.3 a.u. and the initial 
number of protons is 3 141 929. All these spatial distributions contain a 
central maximum and three pairs of symmetric peripheral maxima. In the 
cases of embedded nanotubes the peripheral maxima are labeled by 1^, 2^ and 
3d, and in the case of nanotube in vacuum the peripheral maxima are labeled 
by Ij, 2j and 3j. One can see that the angular distributions obtained with 
the embedded nanotubes do not differ much from each other, and that they 
differ from the angular distribution obtained with the nanotaube in vacuum. 
The spatial distributions obtained with the nanotube embedded in Si02 and 
Ni almost coincide. Figure ITTT d) shows the mappings of the Xq axis in the 
entrance plane of the nanotube to the x axis in its exit plane corresponding 
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to the spatial distributions shown in Figure ITTT c). One can observe that all 
these mappings include three pairs of symmetric extrema. In the cases of 
embedded nanotubes the extrema of the mappings are labeled by 1^, 2^ and 
3d, and in the case of nanotube in vacuum the extrema are labeled by Ij, 2j 
and 3j. In each case the ordinate of the extremum of the mapping coincides 
with the abscissa of the corresponding peripheral maximum of the spatial 
distribution. Therefore, these maxima are the rainbow maxima. 

Let us now compare the above described angular and spatial distributions 
of channeled protons, obtained for v = 5 a.u. and L = 0.5 /im, and for v = 
8 a.u. and L = 0.8 /xm. All of them are characterized by the same proton 
dwell time, i.e., the same duration of the process of proton channeling. When 
the nanotube is in vacuum, the numbers of rainbow maxima of the angular 
and spatial distributions are larger in the former case than in the latter case. 
This is attributed to the weakening of the image force acting on the proton 
with the increase of v. On the other hand, when the nanotube is embedded, 
the numbers of rainbow maxima of the angular and spatial distributions is 
the same in the former and latter cases. This persistency of the rainbow 
effect with the increase of v is explained by the increase of the image force 
due to the polarization of the dielectric boundary. 

Our findings indicate that a dialectric medium surrounding carbon nan- 
otubes can influence their properties. Also, the changing of the spatial distri- 
bution of channeled protons with their kinetic energy we have demonstrated, 
based on the changing of the dynamic polarization effect, may be used to 
probe the atoms or molecules intercalated in the nanotubes. Besides, we 
think that nanotubes could be employed to produce nanosized ion beams for 
applications in biomedicine. 

6.3 Donut effect with protons and a short nanotube in 
vacuum 

In ion channeling experiments the always present questions are the questions 
of ion beam misalignment and divergence. Therefore, it is important to 
explore how the dynamic polarization effect influences the angular and spatial 
distributions of ions channeled through carbon nanotubes when the initial ion 
velocity vector is not parallel to the nanotube axis. We are going to focus on 
the angular distribution of protons channeled through a (11, 9) single-wall 
carbon nanotube in vacuum for v = 3 a.u. and L = 0.2 fim and for the 
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Figure 16: (a) The angular distributions along the G^^ axis of protons trans- 
mitted through a (11, 9) single-wall carbon nanotube embedded in Si02 and 
in vaccum for v = 5 a.u. and L = 0.5 /xm - solid and dashed lines, respec- 
tively, (b) The corresponding deflections functions, (c) The corresponding 
spatial distributions along the x axis in the exit plane of the nanotube - solid 
and dashed lines, respectively, (b) The corresponding mappings of the Xq 
axis in the entrance plane of the nanotube to the x axis in its exit plane - 
solid and dashed lines, respectively. 



angle of the initial prot on velocity vector relative to the nanotuge axis, (p, 
between and 10 mrad (IBorka et al.l . l2010[ ): the corresponding proton kinetic 
energy is 0.223 MeV. In this case the critical angle for channeling is = H 
mrad. It should be noted that the inclusion of a dielectric medium around 
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X (a.u.) Xo (a u ) 



Figure 17: (a) The angular distributions along the axis of protons trans- 
mitted through a (11, 9) single- wall carbon nanotube embedded in Si02, 
AI2O3 and Ni and in vacuum for w = 8 a.u. and L = 0.8 fim - solid, dotted 
and dash dotted lines and dashed line, respectively, (b) The correspond- 
ing Qx{x) deflections functions, (c) The corresponding spatial distributions 
along the x axis in the exit plane of the nanotube - solid and dashed lines, 
respectively, (b) The corresponding mappings of the xq axis in the entrance 
plane of the nanotube to the x axis in its exit plane - solid and dashed lines, 
respectively. 

the nanotube will modify the results of the calculation only slightly. 

It is well-known that when ip is close to the angular distribution 
of ions transmitted through crystal channels takes the shape of a donat 
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(jChaddertonl . Il97nl : iRosner et all . Il978l : lAndersen et all . \l9M). This effect 
is called the donut effect. It has be e n explained afterwards b y the theory 
of crystal rainbows ( Neskovic et al. I . l2002l : iBorka et al.l . 12003 ). It must be 



noted that the donut effect has b een observed in coniputer simula tions of ion 
propagation through nanotubes (IZhevago fc Glebovl . Il998l . l2000h . However, 
the authors did not connect the obtained results to the rainbow effect. 

Figure [TST a) shows the an gular distribution o f transmitted protons for ip 
= 10 mrad, being close to (IBorka et al.l . l2010l ). The image force acting on 
the proton is taken into account. One can see clearly that its has a ring-like 
structure, and that the yield of protons along the ring is very high. This 
is the fully developed donut effect. The corresponding rainbow lines in the 
transmission angle plane are shown in Fig. [THT b). These lines demonstrate 
that the structure of the angular distribution, including the donut effect, is 
due to the rainbow effect. We have established that the influence of the 
dynamic polarization effects on the angular distribution is significant for ip 
between and \l/r7/ 2. However, when g? becomes close to \l/(7, this influence 
is reduced strongly ( Borka et al. . 2010l ). This means that the contribution of 
the dynamic polarization effect to the donut effect is minor. 
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Figure 18: (a) The angular distribution of protons transmitted through a 
(11, 9) single-wall carbon nanotube for v = 3 a.u., L = 0.2 fim and if = 10 
mrad. The image force acting on the proton is taken into account, (b) The 
corresponding rainbow pattern in the transmission angle plane. 
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7 Dynamic polarization effect in proton chan- 
neling through short nanotubes 



The channeling star effect was observed for the first time i n the early exper - 
iments on ion transmission through axial crystal channels ( iGemmelll . Il974j ). 
The characteristic "arms" of a channeling star were attributed to the ions 
moving along the planar channels that cross the axial channel. It must be 
noted that such a motion can occur only whe n the ion bearn . diver gence is 
larger than the critical angle for channeling ( iPetrovic et al.l . |2002| ). Usu- 



ally, the effect has been used for precise alignment of the ion beam with the 
channel axis. 

We are going to describe here the channeling star effect in the transmis- 
sion of protons throu gh the bundles of (10, 10) single-wall carbon nanotubes 
(IBorka et al.l . l2006bl ). The initial proton energy is 1 GeV and the bundle 
length is varied between 2.2 to 14.3 fim, corresponding to the reduced bun- 
dle length associated with the protons moving close to the centers of the 
triangular channels of the bundle, A2, between 0.15 and 1, respectively (see 
part §4 of this review). The study is focused on the problem of mutual 
orientation of the neighboring nanotubes w ithin the bundle, which has not 
been resolved yet ( Kwon fc Tomanek . 2000l ). The nanotubes in question are 
achiral. It is assumed that they form a bundle whose transverse cross-section 
can be descirbed v ia a hexagonal or r hombic superlattice with one nanotube 
per primitive cell ( iThess et al.l . Il996l ). We choose the y axis of the reference 
frame to coincide with the bundle axis and its origin to lie in the entrance 
plane of the bundle. The arrangement of the nanotubes is such that their 
axes intersect the x and y axes of the reference frame, whic h are the verti- 



cal and horizontal axes, respectively (IPetrovic et al.l . l2005d ). We take into 



account the contributions of the nanotubes lying on the two nearest rhom- 
bic coordination lines, relative to the center of the (rhombic) primitive cell 
of the superlattice. The initial proton beam axis coincides with the bundle 
axis. The proton beam divergence angle is chosen to be Qd = Gipc, where 
ipc = 0.314 mrad is the critical angle for channeling. The calculations are 
performed using the theory of crystal rainbows, which is described in part 
§3.2 of this review. The interaction of the proton and a n anotube atom is 
described by the Moliere's approximation (IMolierel . Il947l ) of the Thomas- 
Fermi interaction potential, which is given by Eq. ( !T3|) . The application of 
the continuum approximation gives the Moliere's proton-bundle continuum 
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interaction potential, which is defined by Eq. ( 1T411 . The number of nanotubes 
within the bundle / = 16, the number of atomic strings of a nanotube J = 
40, and the distance between the atoms of an atomic string d = 0.24 nm. 
The thermal vibrations of the nanotube atoms are taken into account. This 
is done by Eq. f|T7|) . The one- di mensional thermal vibration amplitude of 



the nanotube atoms ath = 5.3 pm (IHone et al.l . l2000l ) 



The electronic proton energy loss and dispersion of its proton transmission 
angle, caused by its collisions with the nanotube electrons, are taken into 
account. For the specific proton energy loss we use expression 



2„4 



dEe inZfe 



dz mev'^{z) 



Up 



In 



2mf,'^^v {z) 



(23) 



where me is the electron mass, v{z) the magnitude of the proton velocity 
vector, (3 = v{z)/c, c the speed of light, 7^ = (1 — f3)~^, Ue = {d^x + 
dyy)U^^{x, y)/4:iT the average along the z axis of the density of the nanotube 
electron gas, U^^{x,y) the Moliere's proton-bundle continuum interaction 
potential with the effect of thermal vibrations of the nanotube atoms taken 
into account, h the reduced Planck constant, and = (47re^ne/me) ^ the 
angular frequency of the proton induced oscillations of the nanotube electron 
gas. For the specific change of the dispersion of the proton transmission angle 
we use expression 



dQl 



dz m'iv'^{z) 



dEe 
dz 



(24) 



where irip is the (relativistic) proton mass. The corresponding standard de- 
viations of the components of the proton transmission angle, and Gy, are 



The nanotube diameter is 1.34 nm (jSaito et al.l . 120011 ). and the distance 



between the axes of two neighboring nanotubes is 1.70 nm (IThess et al. 



19961). The diameter of the bundle is chosen to be 169.64 nm, corresponding 
to 100 nanotubes lying on the y axis. The components of the proton impact 
parameter vector are chosen uniformly within the region of the bundle. For 
each proton impact parameter the x and y components of its initial velocity, 
vqx and foy, respectively, are chosen within the Gaussian distributions with 
the standard deviations VLi,x = ^by = ^b/V^, where fib is the proton beam 
divergence angle. The initial number of protons is 859 144. 
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Figure 19: The angular distributions of 1 GeV protons transmitted through 
the bundles of (10, 10) single-wall carbon nanotubes of the reduced lengths 
(a) Aa = 0.15, (b) As = 0.35, (c) As = 0.65, and (d) As = 0.85. The proton 
beam divergence angle is Qd = 1-884 mrad. 

Figure [19] show the angular distributions of transmitted protons for As = 
0.15, 0.35, 0.65 and 0.85. For As = 0.15 a channeling star is clearly visible 
in the angular distribution, and one can distinguish its 40 arms. One can 
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connect this to the fact that each nanotube has 40 atomic strings. This means 
that the channehng star effect can be used to learn about the transverse 
structure of a nanotube. For A2 = 0.35, 12 arms of the channehng star are still 
visible, while for A2 = 0.65 and 0.85, eight arms remain, with the intensity 
being lower in the latter case. The analysis shows that the observed intensity 
lowering and disappearance of the arms as A2 increases can be attributed to 
the process of proton dechanneling, which is very pronounced since Qd ^ ipc- 
For A2 > 1 the channeling star does not exist any more. 

Let us now concentrate on the channeling star for A2 = 0.65. It has two 
arms lying along the 0^; axis, two arms lying along the Qy axis, two arms 
lying along the line defined hy if = tan^^(By/Ba;) = nvr/S and 77r/6, and 
two arms lying along the line defined by ip = 57r/6 and ll7r/6. Figure [20] 
show the components of the proton impact parameter vectors, within the 
(rhombic) primitive cell of the superlattice, that correspond to the arms of 
this channeling star. They are selected taking that the proton transmission 
angle = (B^. + Qy)^ > 1.6 mrad. The analysis shows that these impact 
parameter vectors are located between seven classes of pairs of parallel lines, 
representing seven classes of pairs of planes defined by the atomic strings of 
the nanotubes, which are parallel to the bundle axis. There are two equivalent 
pairs of planes designated by 1, one pair designated by 2, and one pair 
designated by 3, which are shown in Fig. l20r a). and two equivalent pairs of 
planes designated by 4, four equivalent pairs designated by 5, four equivalent 
pairs designated by 6, and four equivalent pairs designated by 7, which are 
shown in Fig. l20r b). Each pair of planes defines a planar channel, which 
corresponds to two opposite arms of the channeling star. The pairs of planes 
that are parallel to each other correspond to the same two opposite arms. In 
accordance with this, the arms shown in Fig. [TW c) are designated by 4, by 
1, 5, 6 and 7, by 2, and by 3. 

In order to analyze the sensitivity of the channeling star effect to the mu- 
tual orientation of the neighboring nanotubes within the bundle, we rotated 
each nanotube about its axis by 7r/60 counterclockwise relative to its assumed 
position 4.2 . The resulting angular distribution of transmitted protons, for 
A2 = 0.65, is given in Fig. [2T1 In this case the channeling star also has eight 
arms, as the one given in Fig. [TW c). but the arrangement of its arms is differ- 
ent. Two of the arms lie along the Qy axis, two arms lie along the line defined 
hj if = n /6 and Ttt/G, two arms lie along the line defined by = 27r/6 and 
= 87r/6, and two arms lie along the line defined by = 57r/6 and lln/G. Fig- 
ure [25] show the components of the proton impact parameter vectors, within 
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Figure 20: The components of the proton impact parameter vectors, within 
the primitive cell of the superlattice, that correspond to the arms of the 
channeling star for the reduced bundle length A2 = 0.65. (a) The solid lines 
represent four pairs of planes defined by the atomic strings of the nanotubes, 
designated by 1, 2 and 3, that enable the planar channeling of the protons 
through the bundle, (b) The solid lines represent 14 additional pairs of planes 
defined by the atomic strings of the nanotubes, designated by 4, 5, 6 and 7 
that enable the planar channeling of the protons through the bundle. 



the (rhombic) primitive cell of the superlattice, that correspond to the arms 
of the channeling star. They were selected taking that > 1.0 mrad. As in 
the above considered case, the impact parameter vectors are located between 
seven classes of pairs of parallel lines, representing seven classes of pairs of 
planes defined by the atomic strings of the nanotubes. There are two equiv- 
alent pairs of planes designated by 1, one pair designated by 2 and one pair 
designated by 3, which are shown in Fig. [2^ a). and two equivalent pairs of 
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A, = 0.65 

h I 




0, [mrad] 

Figure 21: The angular distributions of 1 GeV protons transmitted through 
the bundle of (10, 10) single-wall carbon nanotubes of the reduced length A2 
= 0.65 when the nanotubes are rotated about their axes by 7r/60 counter- 
clockwise relative to their assumed positions. The proton beam divergence 
angle is Qd = 1.884 mrad. 



planes designated by 4, two equivalent pairs designated by 5, two equivalent 
pairs designated by 6, and two equivalent pairs designated by 7, which are 
shown in Fig. l22r b). In accordance with this, the arms shown in Fig. [21] are 
designated by 1, by 2, by 4, and by 3, 5, 6 and 7. 

It should be noted that in both of the above considered cases the pairs 
of planes designated by 1, 2 and 3 do not cross the nanotubes, and that the 
pairs of planes designated by 4, 5, 6 and 7 cross the nanotubes. Namely, 
the protons channeled between planes 1, 2 and 3 move only in between the 
nanotubes while the protons channeled between planes 4, 5, 6 and 7 move 



52 



D. Borka, S. Petrovic and N. Neskovic 




Figure 22: The components of the proton impact parameter vectors, within 
the primitive cell of the superlatticc, that correspond to the arms of the chan- 
neling star for the reduced bundle length A2 = 0.65 when the nanotubes are 
rotated about their axes by 7r/60 counterclockwise relative to their assumed 
positions, (a) The sohd hues represent four pairs of planes defined by the 
atomic strings of the nanotubes, designated by 1, 2 and 3, that enable the 
planar channeling of the protons through the bundle, (b) The solid lines 
represent eight additional pairs of planes defined by the atomic strings of the 
nanotubes, designated by 4, 5, 6 and 7 that enable the planar channeling of 
the protons through the bundle. 



in between as well as through the nanotubes. Also, it is clear that planar 
channels 1, 2 and 3 in the former case practically coincide with planar chan- 
nels 1, 2 and 3 in the latter case. However, planar channels 4, 5, 6 and 7 in 
the former case considerably differ from the planar channels 4, 5, 6 and 7 in 
the latter case. This explains the difference between the arrangements of the 
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arms in these two cases. Thus, one can say that the channehng star effect 
is sensitive to the mutual orientation of the neighboring nanotubes within a 
bundle. This means that it is possible to measure the effect and employ the 
obtained results to determine this orientation. 
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